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ABSTRACT 

In this paper, we introduce the octopole moment of the hght distribution in galaxies as a probe 

of the weak lensing shear field. While traditional ellipticity estimates of the local shear derived from 

" " ■ the quadrupole moment are limited by the width of the intrinsic ellipticity distribution of background 

^— N . galaxies, the dispersion in the intrinsic octopole distribution is expected to be much smaller, implying 

f— ^ ' that the signal from this higher order moment is ultimately limited by measurement noise, and not by 

^Sj . intrinsic scatter. We present the computation of the octopole moment and show that current observations 

I ■ are at the regime where the octopole estimates will soon be able to contribute to the overall accuracy of 

^ ' the estimates of local shear fields. Therefore, the prospects for this estimator from future datasets like 

^~~> , the Advanced Camera for Survey and the Next Generation Space Telescope are very promising. 

^3 ' Subject headings: 

^H \ 

1. INTRODUCTION 

J> ■ The analysis of weak shear fields of background galaxies lensed by clusters has proved to be a rich field of inquiry 

t^^ ; (Blandford & Narayan, 1992; Kaiser, Squires & Broadhurst, 1995, hereafter KSB; Kaiser & Squires, 1993; MeUier 1999; 

OO . Bartelmann & Schneider, 2001 and references therein for applications of weak lensing analysis to observations). However, 

the underlying assumption of these analyses has always been that all of the lensing information can be extracted from 

combinations of the quadrupole moments of the light distribution, which define an ellipticity and an orientation for the 

lensed arclets. Under the assumption that the intrinsic orientation of the lensed galaxy is random, the expectation value 

'^^ ] of the complex ellipticity can be related to the reduced shear. Do combinations of the quadrupoles contain all the useful 

-«^ , information on the local shear field? Refregier & Bacon (2001) point out that any even moment of the light distribution 

i-Q will be preferentially oriented along the shear even in the weak lensing limit. 

Mh^ Moreover, it is generally taken for granted that strong lensing of galaxies, whether producing multiple images or strong 

Q ■ arcs, should be analyzed in an fundamentally different way from weakly lensed "arclets". Kneib et al. (1993), Bartelmann 

?H ■ (1995), and others describe the analysis of a strongly lensed arc which lies near a critical curve, and how such images can 

c/3 , be used to reconstruct the lens potential. 

Ci It is not clear that the transition between arclet and arc is as abrupt as the current mode of analysis would suggest. 

^ Another way of thinking about strong arcs is that lensing induces large octopole moments in the light distribution. This 
octopole has two components: a skewness in the light distribution away from the lens, and a term which expresses itself 
as an arc tangential to the surface of constant shear. Ideally, there is a transitional regime, one in which the shear is 



o 



>< 



JH I sufficiently weak that the image would not be characterized as an arc, but nevertheless, an octopole moment might be 
. . I measurable. 

The potential benefits of constructing this higher order estimator are significant. Despite the anticipated noisiness of the 
octopole signal, an optimal combination of this with the quadrupole signal could significantly improve the overall signal 
to noise, since the absolute value of the intrinsic octopole moments are expected to be small compared to the quadrupole 
moments. Moreover, since the induced octopole moments implicitly probe the radial variation of the shear field, one may 
also get an independent local estimate of the radial derivative of the reduced shear. 

Forthcoming instruments, such as the Next Generation Space Telescope (NGST), ongoing surveys, including the Deep 
Lens Survey (DLS), and current and future cluster observations by the Hubble Space Telescope (HST) using the Advanced 
Camera for Survey, will produce an avalanche of data which need to be fully exploited. The high-quality nature of the 
data means that the limitation may due to large the number of galaxies (crowded fields making shape measurements 
tricky) rather than the quality of the imaging, and thus, the octopole approach to data analysis may be quite fruitful. 

Our approach in this paper is as follows. The basic lensing formalism of linear distortions and the notation used in 
this paper are outlined in § 0. In § H octopole moments are defined, in the context of them arising from the second-order 
effects in a circularly symmetric lens. In § 0, the error analysis is discussed, and it is shown that a combination of the 
quadrupole and octopole can be used to estimate the reduced shear and its spatial derivative. Finally, in § H, we present 
a discussion of future prospects, including application to observational data-sets. 

2. BACKGROUND 

We begin by defining our convention and reminding the reader of the standard weak lensing equations. An excellent 
and exhaustive review of this material can be found in Bartelmann & Schneider (2001), from which our conventions are 
borrowed. 
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2 The Galaxy Octopole Moment as a Probe of Weak Lensing Shear Fields 

A point mass, M, at a distance, Di, will deflect a light beam coming from a source at distance, Dg. In particular, if 
the source is observed at a 2-d angle from the lens, then the beam would be deflected by an angle, a: 

{Di-Ds)Di 4GM^ 
"^^ = -Ds ^^- W 

More generally, if the lens is an extended object, with a surface mass density, S(0), then a dimensionless mass, the 
convergence, may be defined as: 

^^0) ^ M ^ (A - Ds)Di ^nCm _ ^2) 

The convergence may be thought of as a source term for a potential, ip{0)^ and related via a Poisson-like equation: 

s/^m = MO) , (3) 

where all gradients and divergences are calculated in the two-dimensional 0-space. 

Extension of equation (^ to the continuous case, and combination with equation (g) yields a deflection angle, 

a{e) = - I d'e'Kie') ^'^' . U) 

Try >-0'|2 ^ ^ 

A beam observed at angle 6 must therefore have originated in the source plane at 

m = e-a{d), (5) 

under the thin lens approximation. 

Suppose we observe a galaxy with its center of light at the position /3o in the source plane. Throughout this paper, we 
will define a local set of coordinates such that /3o, and correspondingly, 9q (the position of the source in the lens plane) are 
at the origin. This assumed convention does not change the final results, but merely makes the equations more compact. 
Solving the lensing equation to linear order, one finds that both Pq and ^o correspond to the center of mapping (i.e. the 
first order lens mapping from source plane to image plane preserves the position of the center of light). 

In linear analysis both Pq and Oq correspond to the center of light. However, in the higher order analysis, discussed 
below, if the position, /3o, in the source plane is lensed to the foreground, the corresponding lensed position will no longer 
necessarily be the center of light. 

Combining equations (0) & (0), the defiection angle, d{6) can be written explicitly as the gradient of the scalar potential, 
d{9) = y gijj- Since lensing conserves surface brightness, a mapping from foreground to background coordinates is sufficient 
to determine a background brightness map from a foreground one (or vice- versa) provided a full knowledge of the geometry 
of the system (cosmology plus the redshifts of the source and lens) and mass distribution of the lens. Thus, we may expand 
around the origin to determine a deprojection operator on a foreground light distribution, which yields the amplification 
matrix, 

di_(^ _ d^m\ ^ ^ 1 - K - 71 -72 



"^^'^-frV^^-^^V -7. l- + 7r)- («) 



Rigorously speaking, this expression is the first term in a Taylor series expansion of the distortion operator. The term 7 
is a complex shear term, representing the anisotropic part of the distortion, with 7 — |7|e^*'^, and the real and imaginary 
parts being denoted with the subscripts, "1" and "2" respectively, as per convention. Using our locally defined coordinate 
system, we have: 

A = M,e, (7) 

Likewise, assuming the absence of a caustic crossing (1 — k — I7I < 0), this can be inverted uniquely to give a projection 
function, 

0^ - A^r.V, (8) 

In this analysis, we focus on expanding this projection operator to the next higher order to derive the octopole moment 
rather than restricting ourselves to the quadrupole moment alone. In general, researchers have treated weak lensing fields 
in the manner described by KSB or its variants (Bacon, Refregier & Ellis, 2000; van Waerbeke et al. 2001). These 
techniques describe the mapping of source-plane quadrupole light distributions to lens-plane distributions, and thus use 
the observed ellipticity and an assumption of random orientation to invert the shear field. In this work, we aim to 
generalize these transformations to the next higher order. Our notation for the n-th order moments of a galaxy is: 

/ «ymN ^ / dxdyl{x, y){x - x)"(?y - y)'" 
Jdxdyl{x,y) 
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where an overline indicates the mean value, where for the quadrupole n + m — 2. We may redefine the surface brightness, 
I{x, y) to any other function of flux without a loss of generality in the above expression. 

The quadrupole moments are generally translated into a complex ellipticity via one of the relations: 



(a;2) - (y2) + 2i{xy) 



{x^) + (y2) + 2^{x^){y^)-{xyY 
or 



(10) 



^ {x^) - {y^) + 2i{xy) 
^ (x2) + (y2) ^ > 

The expectation value of the complex ellipticities around a circular annulus gives an estimate of the combination of 
parameters: 

E{eg)c^g=-^, (12) 

i — K 

Eixe) ^ 2g (13) 

where g is known as the reduced shear. This expression implies a degeneracy between values of the shear and convergence. 
It should also be noted that these relations hold only in the limit of small variance in the intrinsic ellipticity distribution 
and for small values of g. Strictly speaking, these approximations are invalid in the strong lensing regime characterized 
by the dramatic arcs and multiple images. 

3. THE OCTOPOLE MOMENTS 

In this section, we address the issue of analysis of higher order moments of the light distribution. In particular, we 
show that, even in the weak lensing limit in which one would strictly not expect to be able to detect gravitational arcs, 
an arc-like signature can be found from the octopole moments of the galaxy light distribution. 

In order to simplify this discussion somewhat, we assume a radially symmetric potential throughout, and for a fiducial 
galaxy we further choose the coordinate system with a configuration such that it lies along the positive x-axis. This 
is entirely equivalent to defining a radial and tangential component. However, for practical reasons, we chose to do all 
calculations in Cartesian coordinates. 

Under these assumptions, using only linear theory, the lensed (subscript 6) and unlensed (subscript /3) moments can be 
related in a very straightforward way since the Jacobian is diagonal. Thus, 

{x-y^e = [Ar^r [^^-^r (-"y")/. - (,_J_^)„(i_^'^^).„ (^"^")/^ • (14) 

Application of this transform to the quadrupole terms in the weak lensing limit, yields the complex ellipticity transforma- 
tions above. If galaxies are randomly oriented, as is assumed (however, see Crittenden et al. 2001 for an estimate of the 
degree of expected intrinsic alignments), then the expectation value of any of the intrinsic octopole moments (and hence 
the lensed octopole moments in a circularly symmetric potential) will necessarily vanish. 

However, even in the absence of an intrinsic octopole moment, the second order Taylor expansion of the lensing equation 
can give rise to octopole moments. Consider the limit when the source image has a small but finite size, the lensed field 
may be approximated as: 

(3, ~ A.jOj + -A.j.kOjOk . (15) 

By inspection, the following symmetries: Ai^^j = Aij,k and Aji,k — Aij,k, hold. The corresponding terms (expressible as 
local derivatives of the shear and convergence field) are: 

27 27 

^11,1 = -K.i - 71.1 = -271,1 - 72,2 = -271,1 = -27' (16) 

r r 

^12,1 = ^21,1 = ^11,2 = -72,1 = 

27 

^22,1 = ^12.2 = ^21,2 = -K.i + 7i,i = —72.2 = 

r 

-422,2 = -K,2 +71.2 = , 

where 7' — d^/dr, the radial derivative of the shear field. 

Since the potential, and thus the surface density of the lens, S(6'), is circularly symmetric, and the center of the lens 
lies on the x-axis, k,2 — 0. Likewise, since everywhere on the x-axis, 72 — 0, 72,1 = 0. In addition, we have used the 
relation derived by Kaiser (1995): 

( '^A \ - f Tia +72,2 
Finally, it can be shown that 72,2 — 27/r. 



K,2 J V 72,1 - 7l,2 ' ^ ' 
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Summarizing and incorporating all the above symmetries and simplifications, we have: 



7 



2 7 2 

r 



X0 ~ Aiixe - 7' + 

L /■ 

A 27 
yp ^ ^222/9 Xeyg 



Thus, to second order in the position, the above equations may be inverted: 



(18) 



xe ~ Aj^^xp + Aj^{^ 



7' + 



7 



+ ^11 ^22 ^y/3 



(19) 



ye - A^^'yp + A^^'A 



-1.-227. 



22 — -^py^ 



In order to compute the moments of the lensed field, we need to transform the area element by determining the Jacobian 
consistently to the same order. 



J; 



dx0 dxg 
9xf3 dyf, 
dxg dxg 
dx/j dyf, 








^n + 2^n [7' + ?] xp 2A-lA-^:Ly^ 

2ArM2t ^2//3 A22 + 2An ^2^ ?^/3 


u + {2A:[^A22 




+ 2Ari%tJ)2^;3, 



(20) 



to first order in Xfj. This yields an area clement in the lensed field such that dxgdyg =^ Jdxpdyp. 

We may now compute the expected transformation of any given moment. The total flux, / from the galaxy, for example 
is: 



fe = / I{x,y)Jdxpdyp = ^/^ 



(21) 

which is precisely the result in linear theory. We may also compute the shift in the center of light compared to the 
linear theory prediction. Symmetry arguments, along with our choice of coordinates clearly will produce no shift in the 
y-coordinate direction. 
However, 



{x)e 



1 



M//3 



[a^^xp + A^^ [y + J] xl + An'^22'^y') Jdxpdyp 



(22) 



2A-,^ 7' + - + '^A-^iA^^l {x^)p + A^lA2^^{y 



2 7, ,2\ 



l-lTl 



-l7/„.2\ 



3^n (7' + ^)+ 2A2ij:\ {x')e + A2iMy^)e 

Since this shift merely reflects a change in position between the flrst and second order approximation of the lensing 
inversion, it is not directly measurable. On the other hand, we may compute the expectation values of the observables, 
{x^)0 and {xy'^)g: 



{x')e = Ari'(^'>/3 + UA22'- + 5A^i (7' + -)1 {x*)e + 3A^,'l{x^y\ 
9An (7' + ^)+ eA^fA {xYe iA-,l^-{x^)o{yYe 



(23) 



-3/„3\ 



AV'lx' 



r(l + g) 



9' - g/r - g^/r 



6 



.9' ^ 9/1" ^ .9^/^ 
I-.92 

2X2 .3 



{x^)e + 3 



r(l + g) 



{x~ 



Ki - .9) 



r{l - .9) 

x'MyYe 



{xhh 



and 



{xy^)g = A^^A2^{xy^)p + [dA'^^ (7' + J) + QA^^-] {x'y\ + ^n J(2/'>e 
- Mri^ (7' + ^)+ 2A-f-\ {x\{y')g A-f-{yYe 



(24) 



AiiA22{xy'^)p + 

9' -glr-g^l 
1-52 



9' - g/r-g'^h 



l-.g2 



(-^y^>. 



1 + 5 



{y')B 



1 + 5, 



(^'>«(5'> 



1-5 



{y 



2\2 



Goldberg & Natarajan 5 

In both these expressions (eqns. 23 and 24), the terms proportional to squares of the quadrupole moments arise from a 
second order shift in the center of hght. By inspection, the expectation values of the other two octopole moments, (x^y)^ 
& (y'^)e, will vanish. 

Note that with the exception of the term proportional to the intrinsic octopole moments, both the lensed octopoles are 
a function of observables: the second and fourth moments of the light, the reduced shear, gr, and its first derivative. The 
octopole analysis, however, retains the degeneracy found in quadrupole analysis. However, since the scatter in the intrinsic 
octopoles may be much smaller than the scatter in the intrinsic ellipticity, and since it has a mean of zero, extracting a 
profile of g(r) from a finite number of noisy galaxies may be done with potentially higher signal to noise by using the 
third as well as the second moments. We discuss this parameter extraction in the next section. 

4. ERROR ANALYSIS AND PARAMETER ESTIMATION 
4.1. Estimation of the Errors on the Moments 

The fact that lensing induces an octopole moment in an otherwise elliptical galaxy is interesting, but by no means 
useful unless we are actually able to measure this effect. We must measure both the octopole moments and the hexade- 
cipole moments with some accuracy in order to extract the shear, and since both moments are much more likely to be 
contaminated by sky noise than their quadrupole counterparts, it is possible that measuring these moments may be quite 
difficult. In order to estimate the feasibility of such an investigation, we must first consider the uncertainties in measuring 
the moments for a fiducial galaxy. 

We make a number of simplifying assumptions in estimating the moment uncertainties. First, we assume that the only 
source of error arises from the Poisson noise in the galaxy image and the background sky. We thus do not include the 
effects of a PSF, which are likely to be considerably more important on the lower order quadrupole moments than on the 
octo- and hexadecipoles. Similarly, we do not include pixelation effects, which are also expected to be more important 
for the quadrupolcs. However, we do assume that the background sky level is constant (up to Poisson noise) over the 
entire image, and that the galaxy image is unblended with any of its neighbors. It is these effects that will contribute 
significantly to the noise in the higher order moments. 

Given these constraints, we assume that a good estimate of the moments of a galaxy can be achieved by integrating using 
a circular mask, including a tophat and Gaussian. Although, this is clearly not optimal, as a circular mask will introduce 
biases in all moments, we adopt this aperture for purposes of illustration in this work. The generality of the formalism 
developed here allows one to potentially replace the surface brightness in equation (0), with any other monotonic function 
of surface brightness. However, optimal determination of this function would rely on some knowledge of the true radial 
light profile of the galaxy. 

Clearly, the analysis can be successfully performed for an elliptical mask, or even an arbitrarily shaped mask, in which 
instance the relevant region above some given isophote could be determined iteratively. Discussion and application of 
optimal shape estimation can be found in the literature on IMCAT (Kaiser), sExtractor (Bertin & Arnouts, 1996), FOG AS 
(Jarvis & Tyson 1981). However, these pieces of software have focused on estimating only the ellipticity and other shape 
parameters related to the quadrupole moments. 

Specially promising among recent shape estimation methods is the Shapelets technique (Refregier & Bacon, 2001; 
Refregier, 2001) which decomposes images into Hermite polynomial basis sets. A potential approach to computing 
octopole moments may be to determine the second order shear operator directly in shapelet space. 

However, it is beyond the scope of the analysis of this work to investigate optimal techniques for measuring higher order 
moments. This will be dealt with in future work. For the present, we will simply use the assumptions above to effectively 
characterize the practical measurement uncertainties in typical images. 

That said, one may approximate the uncertainty in (a;"j/™), an,m as: 

2 _ /o " de cos2"(0) sin2"'(0) J^ dr w^{r, 6) ri+2«+2™ [/(r, 9) + N] 



{Cd(^Io°°drw{r,e)rI{r,e) 



where /(r, 9) is the surface brightness in counts/area of the galaxy, N is the background sky brightness in counts/area, 
and w{r, 6) is a mask. For our analysis, we use two forms for this mask, a tophat, with w = 1 out to a specified outer 
radius, and a circular Gaussian, as used in KSB. 

The form of the integral explicitly assumes that the light is very nearly circularly distributed. For an assumed radial 
profile, the error may be approximated analytically or semi-analytically. For our analysis, we have used a de Vaucouleurs 
(1948) profile: 

/(r, 6) = I, exp{-7.67[(i?/i?,)i/4 - 1]} , (26) 

where le is the characteristic brightness, and Re is the half light radius. We find similar error estimates for other assumed 
profiles. Typically, we set the outer radius at 4 times the half-light radius. 

Note that we ignore any covariances in the errors. While inclusion of these terms in our analysis would generally lower 
the implied errors, their effect is expected to be small compared to other uncertainties which we have explicitly overlooked. 

In addition to measurement errors, both the ellipticity estimate of the shear and the octopole measurement of the shear 
require an estimate of the intrinsic scatter of the ellipticity or the octopole. The first we will label cr^, which is somewhat 
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different from the value normally given. Since we assume a random orientation, a^ represents the standard deviation in 
the intrinsic distribution of the real part of the complex ellipticity, a term which has an expected mean of zero. 

The intrinsic variation in the octopole moments in the octopoles are labeled cr(2:3), and ai^^yiy These terms represent 
the standard deviation in the following dimcnsionless form: 

/ 3\\2\ ^/^ 

(27) 

with a similar form for the other octopole. We again assume no covariance. 

To simplify the analysis, we assume that the half-light radius, i?e, and the characteristic radius is known with perfect 
certainty, as is the integrated flux. Since these terms are dominant near the center of the image, their errors will be 
considerably smaller than the higher order moments. 

4.2. Uncertainty from Ellipticity Estimates 

Finally, we consider the uncertainty in the shear from the measurement of the quadrupoles of a single galaxy. In the 
limit of weak lensing, the reduced shear may be approximated as: 

_ Xi - Xi'' ^ jx^) - jy^) _ Xi^ .28) 

^ 2 2r2 2 ■ ^ ^ 

Since the orientation of the intrinsic ellipticity, x^ is random, ct^^^ = CT;^/v2- Ebbels et al. (2000) derive a probability 
distribution function for the observed ellipticity distribution for galaxies of different morphological types from the Medium 
Deep Survey. Adapting this, we find a reasonable value of a^ = 0.30. 
Thus, the variance in the estimate of g from the ellipticity is: 

In most instances, this uncertainty will be dominated by the spread in the intrinsic galaxy ellipticities. 

4.3. Uncertainty from Octopole Estimates 

The details of parameter estimation from the octopole moments are slightly more involved. We may begin by simplifying 
the terms in eqns. (p4 and Eq) which are functions of g. For small g, we may say: 

9 rg' -g-g^ , 

TT— -^ - 9 ; ] 2 — -^9-9 (30 

(1 + 5) 1-5^ 

We then simplify eqns. (|2J and |2^) to the following expressions: 

(a;^)e = ^n (a;^>/3 + Uii5 + "127- 9' , (31) 

and 

{xy^)e = A];^A22{xy^)p + '"2i.g + ^22?^ g' , (32) 

where the terms, Uij, can be derived by inspection, and their uncertainties can be determined by the standard propagation 
of errors. We reiterate here that this analysis assumes that the errors are Gaussian and uncorrelated, which is certainly 
not the case. 

Inverting these expressions yields: 

9 [u22i{xy^)e - A^^^ A^Hxy^) ,3) - ui2((V>e - A^^^ A^^ (xy^) p)] (33) 

U11U22 — U12U21 

and 

r 9'- [-U2i{{xy^)e - A^^^ A22^ (xy^) p) + uiii{xy^)g - A^^ A22^ (xy^) p)] (34) 

U11U22 -U12U21 

Since the expectation value of the intrinsic octopoles is zero, estimating the parameters from these equations is straight- 
forward. 

Finally, from the form above, calculation of the variance and covariance of g and r g' is a straightforward but tedious 
exercise in propagation of errors. 
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4.4. Results 

It now remains for us to compare errors estimated from the quadrupole moments alone to those from the octopole 
moments. For these estimates, we assume the galaxy to intrinsically follow a de Vaucouleurs profile. Note that this 
assumption is necessary only in computing the estimate of the error, and is not required for parameter estimation. 

The fiducial galaxy has a magnitude of 22 at 2; = 0.5, and has a half-light radius in the source plane of 2kpc. Error 
estimates are made for observations taken using the HST Wide Field camera using a 16.8 ks exposure as was used for 
observations of AC114 (Natarajan et al. 1998). Signal and sky counts were estimated using the HST WFPC2 exposure 
time calculator []. The lens and source are separated by 100 pixels in the image, and we take a fiducial shear of 2% at the 
center of the source. The lens is assumed to be isothermal, with 2% mean shear at that radius. The precise amplitude of 
the shear is not that important in comparing the ellipticity and octopole estimate error bars, since both techniques are, 
in fact, first order in the shear. 

The deviation in the intrinsic ellipticity, a^ is taken as 0.30, as discussed above while the scatter in both the intrinsic 
octopoles is 0.1. The intrinsic octopole distribution is not well studied, and understanding this may prove to be a fruitful 
inquiry into galaxy shapes. We select a very small scatter here for the following reason. If the intrinsic dispersion in 
the octopole moments is sufficiently large that for most observational sets the uncertainties in parameter estimation are 
dominated not by photon Poisson noise, but by shot noise in the sample, then the octopole moment at best improves the 
signal to noise compared to the quadrupole by v2- The transition from Poisson to shot noise domination is illustrated in 
Fig. 2g. 

The error ellipses from the quadrupole and octopole estimates of g and rg' , along with their combined ellipse can be 
found in Figure |l|. Note that even for existing data sets, inclusion of the octopole analysis has the potential to yield 
significant new information about the shear field. 

The effect of the choice of physical and observational parameters on the estimated errors are shown in Figure |l|. The 
error ellipses for a slice through parameter space are shown in Figure g. Analysis is performed using a Gaussian mask 
with an effective radius equal to the half light radius. 

Note that variations in the observational and physical parameters within these ranges do not generally alter the error 
from the quadrupole moment alone. This is because in this regime, the quadrupole error is dominated by the intrinsic 
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Fig. 1. — The calculated error ellipses using the method developed in this paper. The dotted lines represent the Icr errors on g from the 
quadrupole technique alone. The dashed line represents the 1 — a error ellipse on g and rg' from the octopole method. The solid ellipse is 
the combined error estimate. In each panel, the simulated galaxy has been assigned an apparent magnitude of 22 at 2 = 0.5. Calculations 
are done for the HST Wide Field Camera using the 814nm filter. Exposure time is set for 16.8ks, a^ = 0.30, and the relative variances for 
both relevant octopole moments are set to 0.1. The reduced shear is set to 2%, and the source appears 100 pixels from the lens. The four 
panels use different masks to do the signal to noise estimate. Panels a) and c) use Gaussian masks, as described in KSB, with characteristic 
radii set to the Re and 2/?e, respectively. Panels b) and d) use tophat masks, with radii equal to 2Re and ARe, respectively. 



The Galaxy Octopole Moment as a Probe of Weak Lensing Shear Fields 




0.4 
0.3 
0.2 
0.1 



21 21.5 22 22.5 



e) R (pixels) 



0.2 
0,15 



0.25 


U 






0.2 


r 


b) Source Rcdshift 


- 


0.1b 


— ^ 


-^ 


— 




rl 


1 1 1 1 1 1 1 T~l^i 1 1 1 1 in r 


T—r I- 



0.3 0.4 0.5 0.6 0.7 




.L. J. . L.J. .L.J. . I.. J. 

f) Local shear, y 



50 



100 



150 



0.01 0.02 0.03 0.04 




Quad. Est. of g 
Oct. Est. of g 
Oct. Est. of rg' 



0.05 0.1 0.15 



Fig. 2. — A sequence of slices through observational and physical parameters, comparing the variation of uncertainty in parameter estimates. 
We vary a) exposure time, b) source redshift, c) source apparent magnitude, d) source effective radius, e) source-lens separation on the chip, 
f) local shear on the source, and g) intrinsic scatter in the octopole. Note that the uncertainty in g arising from the quadrupole estimate 
alone remains essentially unchanged over this set of parameters. This is due to the fact that within this regime the dominant effect in the 
quadrupole signal is produced by shot noise. 
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variation in ellipticities, rather than measurement errors. If such is the case, then we are necessarily limited by shot noise 
in the number of sources rather than Poisson noise. If the intrinsic scatter in the octopole moment is sufficiently large, 
this may be the case for the octopole estimator as well. However, we reiterate that this statistic is not well-studied. 

For brighter sources or longer exposures, the quality of estimates for the octopoles improves quickly, as do the estimates 
if the distance to the lens decreases, even if this is not accompanied by a corresponding increase in local shear. 

5. DISCUSSION AND FUTURE PROSPECTS 

In this paper, we have introduced a new and potentially powerful way of analyzing weak lensing shear fields - using the 
octopole moments of the observed light distribution as a second-order estimator of the shear. In addition to the shear, we 
also obtain an estimate of its radial derivative. The radial derivative of the shear field provides an important constraint on 
the mass profiles of the dark matter halos that host galaxies. There are few other methods that can probe that variation. 

While we have demonstrated that within the range of reasonable physical and observational parameters the correspond- 
ing measurement uncertainties for the octopole may be comparable to those found using traditional ellipticity estimates 
of the shear, much remains to be done, both theoretically and observationally. The Gaussian mask used in estimating the 
octopole and quadrupole moments is almost certainly not optimal for this technique. In order to apply this to observa- 
tions, shape estimators must be used which can compute these higher-order moments with maximum signal to noise. One 
avenue of inquiry is to apply a second-order analysis to the Shapelet (Refregier, 2001; Refregier & Bacon, 2001) technique 
in order to get a comparable signal. 

In addition, we have created a somewhat simplistic model of the lens as a circular system. Schneider & Bartelmann 
(1997), for example, consider an extension of KSB with an elliptical lens. Such an effect would introduce an additional 
degree of freedom in the present analysis, since the shear and its radial derivative would no longer necessarily align. 
Application of this technique with generality must include these effects. 

Finally, this technique needs to be applied to data. A number of excellent observational datasets exist with comparable 
observational parameters to those used in the text (e.g. AC114, Natarajan et al. 1998). In addition, the Deep Lens 
Survey |^ (DLS) , and the forthcoming HST Advanced Camera for Survey (ACS) will provide the sort of high-quality data 
which will be ideal for this analysis. We hope to present an application to currently available data-sets in a future paper. 
It is clear that there is significant information in higher-order moments of the light distribution of galaxies, beyond just 
the ellipticity, and here we have illustrated the uses of the octopole moment. 

We would like to thank Ue-Li Pen and Tereasa Brainerd for helpful discussions. 
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